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Abstract 

We calculate the cross sections for the direct hadroproduction of J/tp and T asso- 
ciated with a heavy-quark pair of the same flavour at leading order in as and v 
in NRQCD. These processes provide an interesting signature that could be studied 
at the Tevatron and the LHC and also constitute a gauge-invariant subset of the 
NLO corrections to the inclusive hadroproduction of J/ip and T. We find that the 
fragmentation approximation commonly used to evaluate the contribution of these 
processes to the inclusive quarkonium production sizeably underestimates the exact 
calculation in the kinematical region accessible at the Tevatron. Both J/ip and T 
are predicted to be unpolarised, independently of their transverse momentum. 
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1 Introduction 



Since its introduction in 1995 [1,2], NRQCD has become the standard frame- 
work to study heavy-quarkonium physics. NRQCD is a non-relativistic effec- 
tive theory equivalent to QCD, where inclusive cross sections and decay rates 
can be factorised as the product of short- distance and long-distance parts. 
The short-distance coefficients are process-dependent but can be computed 
in perturbative QCD, the strong coupling constant as being the parameter 
of the expansion. The long-distance matrix elements are process-independent, 
i.e., universal, but non-perturbative. They can be classified in terms of their 
scaling in v, the relative velocity of the heavy quarks in the bound state, and 
as a result, physical quantities can be expressed as a double series in a s and 
v. 
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An innovation of this formalism is the colour-octet mechanism: the heavy- 
quark pair is allowed to be created in a colour-octet state over short distances, 
the color being neutralised over long distances. It is thanks to this very mech- 
anism that it is possible to account for the CDF data [3,4] on the inclusive 
J/ip and ip' cross sections at the Tevatron [5,6,7]. In this case the colour-octet 
matrix elements fitted from the data roughly scale as the power counting rules 
of NRQCD predict. On the other hand, the recent data collected at y/s = 1.96 
TeV by the CDF collaboration [8] have revealed that the J/ip is unpolarised, 
in flagrant disagreement with the expectations of NRQCD. It is fair to say that 
the mechanisms responsible for the quarkonium production at the Tevatron 
are not completely understood yet (for a recent review see [9]). 

Another challenge to theorists has been provided by the recent measurements 
at B factories, where rates for inclusive and exclusive J/ip production in asso- 
ciation with a cc quark pair [10] are far larger than the leading order NRQCD 
theoretical expectations [11]. In this case, keeping only the leading term in as 
and v in the NRQCD expansion results in a very rough approximation. In addi- 
tion, as opposed to the production of the J / ip in hadron collisions, colour-octet 
channels cannot be invoked to bring theoretical predictions in agreement with 
the experimental measurement [12]. Inclusion of the a| (NLO) corrections to 
this process has been recently found to reduce the discrepancy between theory 
and measurements [13]. 

In view of the unexpectedly large measurements for the production of J/ip 
associated with a cc quark pair in e + e~ annihilation, it is natural to wonder 
whether the corresponding production pattern could be large in hadropro- 
duction as well. Besides offering a new interesting signature, such as £~£ + in 
association with one or two heavy quark tags, this process contributes to the 
a s (NLO) corrections to the inclusive colour-singlet hadroproduction of J/ip 
and T. Historically, these higher-order contributions to the cross section at 
the Tevatron were considered in the fragmentation approximation as a first 
attempt to solve the so-called ip' anomaly [14,15]. 

The purpose of this work is to present and discuss the complete tree-level 
calculation for the associated production of J/ip and T (commonly noted 
Q) with a heavy-quark pair of the same flavour. Such associated production 
has been recently discussed in the kx factorisation formalism [16,17]. Here 
we consider the leading order term in the NRQCD expansion in the usual 
collinear factorisation scheme. 

This Letter is organised as follows. In Section 2, we briefly decribe the method 
used to compute the associated hadroproduction of a quarkonium with a 
heavy-quark pair together with the results for the Tevatron and the LHC. In 
Section 3, we discuss the validity of the fragmentation approximation Q — > QQ 
at the Tevatron. In Section 4 we discuss the polarisation of the quarkonium 



2 




Fig. 1. Representative diagrams contributing at LO in as to gg —* Q + via a 
3 S[ 1] state. 

produced by colour-singlet transitions. Finally, Section 5 is devoted to the 
discussion of our results. 



2 The calculation of pp — > Q + 



S.i T/ie method 



In this section, we present the numerical method used to compute quarkonium 
production amplitudes. At leading order in v, the invariant amplitude to create 
a 3 Si quarkonium Q of momentum P and polarisation A can be expressed as 
the product of the amplitude to create the corresponding heavy-quark pair, a 
spin projector iV(A|si, and R(0), the radial wave function at the origin in 
the configuration spacq_j, namely, [18,19] 



M{gg -> Q t Qj(p = 0; s u s 2 ) + X) , 

where P = pq +pq, p = PQ 2 Q , si is the quark spin, s 2 the antiquark one and 
' is the projector onto a color-singlet state. In the non-relativistic limit, the 



1 Up to v 4 corrections, the relation to the NRQCD production matrix element is 
given by {0\)j/^ = 4^|-R(0)| 2 , where the normalisation of Ref. [1] is assumed. 
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spin projector factor can be written as 



N(X\si, to) = ^^«(|. *0 VV§, *i) , (2) 
where e x is the polarisation vector of the quarkonium. 

The amplitude for gg — > Q X (P) + QQ involves 42 Feynman diagrams and an 
analytical computation is not practical. We have therefore opted for using a 
custom version of MadGraph [20] , to generate the amplitudes and perform the 
spin projection numerically. The expression in Eq. (2) can be evaluated with 
the help of the HELAS subroutines [21], as done in Ref. [11]. We validated our 
algorithm to generate the amplitudes for quarkonium production and decays 
by comparing with several known analytical results point-by-point in phase 
space. Finally, we use the techniques introduced in Ref. [22] to perform the 
phase-space integration. As an exercise, we have reproduced the cross sec- 
tion and the P T distribution for B* production at the Tevatron and found 
agreement with the results of Ref. [23] . 

At hadron colliders a 3 Si quarkonium state is observed through its decay into 
leptons. In order to keep the spin correlations between the vector-like bound 
state and the leptons, we can replace e x in Eq. (2) by the leptonic current 

u e -(k 1: Ai)7 M ^+(A;2, A 2 ) . (3) 



8m QV /7r 

The polarisation of the quarkonium can therefore be determined by analysing 
the angular distribution of the leptons. Defining 6 as the angle between the £ + 
direction in the quarkonium rest frame and the quarkonium direction in the 
laboratory frame, the normalised angular distribution J(cos(6>)) is 

/(costf) = (1 + a cos 2 9) , (4) 

where the relation between a and the polarisation state of the quarkonium is 

a T - 2a L 

a = — . (5) 



2.2 Cross sections and results 



As for the case of open charm and bottom cross sections, heavy-quarkonium 
hadropro duct ion is dominated by gluon fusion. We have checked that the light- 
quark initiated process for Q + QQ production is suppressed by three orders 
of magnitude, and thus this contribution is neglected in the following. 
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In our numerical studies we have used: 



• \Rj/40)\ 2 = 0.81 GeV 3 and \R T{1S) (0)\ 2 = 6.48 GeV 3 ; 

• n Q = v /(4m Q ) 2 + p2 ; 

• Br(J/ip -> fj+fi-) = 0.0588 and Br(T(lS) -> //+//") = 0.0248 

• m c = 1.5 GeV and m fe = 4.75 GeV; 

• pdf set: CTEQ6M [24]. 



In Fig. 2 we show the P T distributions of J/ip and T at the Tevatron and the 
LHC, for both Q + g and Q + QQ production. The theoretical error bands in 
the Q + QQ curves correspond to the uncertainties from the renormalisation 
and factorisation scale ( ^ < Hf r < 2/i ) and the heavy-quark masses (m c = 
1.5 ± 0.1 GeV or m& = 4.75 ± 0.05 GeV), combined in quadrature. 



We start by considering J/ip and T production at the Tevatron, where a 
rapidity cut \y\ < 0.6 is applied. For the charmonium case, we note that the Pt 
distribution peaks at Pt — itiq and then it starts a quick descrease, dropping 
by four orders of magnitude at Pt — 20 GeV. In fact, at moderate transverse 
momentum, Pt > mg, Q + QQ production has much milder slope compared 
to that of Q + g. Even though suppressed by as, Q + QQ production starts to 
be comparable to Q + g already for P T — 4:mQ. We verified that the topologies 
where the Q is produced by two different quark lines always dominate. This 
feature is to be related with the J/ip production at e + e~ colliders, where the 
presence of two extra charmed mesons seems to indicate that the J/ is mostly 
often created from two different quark lines as well. For the T, the curve is 
less steep in the same Pt range due to the larger value of the 6-quark mass. 



The results for the LHC are displayed in Fig. 2 (c,d). A cut \y\ < 0.5 has been 
applied. Since the rapidity distribution is flat in this range, this is equivalent 
to consider the quarkonium production at zero rapidity. The Pt behaviours 
for both J/ip and T are very similar to those at the Tevatron while the nor- 
malisation is increased by one order of magnitude. 
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Fig. 2. (a) Differential cross section for the process pp — ► J/ip + cc at the Teva- 
tron, y/s = 1.96 TeV. For comparison, the Pt shape associated to the process 
pp — ► J/tp + g is also displayed (dotted line), (b) Same plot for T + bb production. 
(c,d) Same plot as (a,b) at the LHC, y/s = 14 TeV. 



3 Testing the fragmentation approximation 



In this section we discuss the range of applicability of the fragmentation ap- 
proximation for the process where a QQ pair is produced at high Pt and one 
of the heavy quarks fragments into a quarkonium state, Q — > QQ. At suf- 
ficiently large Pt, it is expected that the quark- fragmentation contributions 
would dominate the QQQ cross section. 



6 



3.1 Heavy-quark fragmentation 



In the fragmentation approximation, the cross section for the production of a 
quarkonium Q by gluon fusion via a heavy quark Qi fragmentation is given, 
at all orders in as, by 

r 1 P 

d<TQ(P)=22 dz da Qr( — t 'rag) Dq.-+q(z, fi f rag ), (6) 
i J0 Z 

where Axq^-j, /i/ m9 ) is the differential cross section to produce an on-shell 
heavy quark Qi with momentum — and Dq.^q(z, Hf rag ) is the fragmentation 
function of Qi into a quarkonium Q. 

The fragmentation scale, fif rag , is usually chosen to avoid large logarithms of 
Pr/n-frag in CTQ^^/ifrag), that is Hfrag — Pt- The summation of the corre- 
sponding large logarithms of Hfrag/mQ appearing in the fragmentation func- 
tion can be obtained via an evolution equation [25,26,27]. 

The perturbative quark fragmentation function into quarkonium Q via a 
state at the scale 3mQ is [28] 

. 8a|(2m Q ) \R(0)\ 2 z(l -z) 2 (16- 32z + 12z 2 - 32z 3 + 5z*) 
D Q ^ Q (z,3m Q ) = -g ^—- 6 . 

(7) 



3.2 Comparison 



We now compare the results of the full LO cross sections for pp — > Q + QQ 
with those calculated in the fragmentation approximation. The same set of 
parameters of Section 2.2 is employed"^!. The comparison is shown in Fig. 3. 
The approximation is lower than the full computation in the Pt range accessi- 
ble at the Tevatron. We have verified that for J/ip production, the two curves 
still differ by a little bit less than 10% at P T = 80 GeV. For the T, much 
larger values of Pt are required (Fig. 3 (b)). 

Therefore, contrary to the common wisdom, in the kinematical region accessi- 
ble at the Tevatron, pp — > J/ip + cc, which is an NLO subset of pp —>■ J/ip + X, 
is not dominated by the fragmentation contributions. 



2 For consistency, the coupling constant in the fragmentation function is evaluated 
at the scale yj (4mQ) 2 + P^,. Note that this choice, as well as the use of CTEQ6M, 
leads to smaller cross sections compared those of Ref . [28] . 



7 





1e-05 



1e-10 



P T (GeV) 

(a) 



20 30 40 50 
P T (GeV) 

(b) 



60 



Fig. 3. (a) Comparison between the full LO cross section for pp — ► J/ip + cc and the 
fragmentation approximation at y/s = 1.96 TeV. No cut on rapidity is applied, (b) 
Idem for pp — > T + bb. 

4 Polarisation: gg — > Q + <3Q vs. — > Q + g 



Let us first review the analytical result for gg — > Q + g, which we used as a 
check of our numerical procedure. To compute the polarisation parameter a 
discussed in Section 2.1, it is sufficient to consider the unpolarised (or total) 
cross section and the longitudinal one in the laboratory frame. The parameter 
a introduced in Section 2.1 can then be written as 



a = 



a T — 2a L _ o to t — 3o"l 
a T J r 2a L cr tot + a L 



(8) 



The total cross section is calculated from the squared amplitude summed over 
the quarkonium polarisation which reads [18,29] 



\M 2 = 



320g 6 s \R(0)\ 2 s 2 M Q f(t + u) 2 + (s + uf + (s + tf 



9tt 



{s + t) 2 {s + u) 2 {t + u) 2 



(9) 



where Mq = 2m<g and s, t and u are the usual Mandelstam variables for 
parton quantities. To compute the squared amplitude for longitudinally po- 
larised quarkonia, we need to define the longitudinal polarisation vector in the 
laboratory frame: 



e\ (P) = a L k x + b L k 2 + c L k A , 



(10) 
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with 



a L 
b L 
cl 



Mq(xi — x 2 ) — (tx 2 + UXi) 
M Q N 

Mq(x 2 - xi) - (ix 2 + uxi) 

M Q N (11) 
-Mg(xi + x 2 ) + (ix 2 + uxi) 
M Q N 

N= J((xi - x 2 ) 2 M^ - 2{u Xl - tx 2 )( Xl - x 2 )M 2 Q + (tx 2 + u Xl ) 2 ) , 

where X\ and x 2 are the usual momentum fractions of the incoming gluons. 
From that definition and the Eqs. (10.505) of Ref. [30], we can then easily 
compute 

QA0g 6 s \R(0)\ 2 M^ Q siu (gsj + s\ Xl - x 2 f + u 2 x\) 
1 Ll 97r(t + M) 2 (s + M) 2 (s + t) 2 iV 2 ' 

If Xi and x 2 are approximately equal and the masses can be neglected, it is 
possible to see from Eq. (9) and Eq. (12) that the quarkonia will be transver- 
sally polarised. Using the same parameters as for the previous plots, we find 
that a gets rapidly close to one for P T > rriQ, see Fig. 4. Our results agree 
with those of Ref. [31]. 




4 8 12 16 20 4 8 12 16 20 

P T (GeV) P T (GeV) 

(a) (b) 

Fig. 4. Polarisation parameter for pp — > Q+g for J/ip (a) and T (b) at the Tevatron, 
y/s = 1.96 TeV. 



Note that this behaviour of a is not expected to be seen in the data, since the 
yield from pp — > Q + g is not dominant, even when only colour-singlet channels 
are considered. For the process gg — > Q + QQ, a remains close to zero in the 
whole P T range. Since the latter is the one which dominates over the other 
colour-singlet contributions, we conclude that the quarkonia not produced by 
a colour-octet mechanism will be unpolarised at high P T . 
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Fig. 5. Polarisation parameter for pp 
Tevatron, ^/i = 1.96 TeV. 
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Q + QQ for J /iff (a) and T (b) at the 



5 Conclusion 



In this Letter, we have presented the tree-level calculation for the associated 
production of J / ip and T with a heavy-quark pair. The motivation was twofold. 
First, these processes offer a new interesting signature, that could be tested 
experimentally by measuring the fraction of quarkonium produced with at 
least one heavy-light quark meson. Second, they contribute to the a| (NLO) 
corrections to the inclusive colour-singlet hadropro duct ion of J ftp and T. 

We have found that the quark-fragmentation approximation employed to de- 
scribe quarkonium production should not be applied in the range of transverse 
momenta reached at the Tevatron and analysed by the CDF collaboration. The 
fragmentation approximation actually underestimates the full contribution of 
pp — > J/ip + cc by more than a factor four in the region P T ~ 10 GeV. This 
approximation starts to be justified (within 10% of error) only for Pt > 80 
GeV. The same conclusion applies for the ip'. For pp — > T + bb the fragmenta- 
tion approximation is relevant at even higher Pt, as expected. It is interesting 
to note that for J/ip production from e + e~ annihilation, the fragmentation 
approximation is only 6% away from the full contribution at yfs = 50 GeV 
[11]. One may try to gather an intuitive explanation for such glaring difference 
by considering the Feynman diagrams contributing to the two processes. In 
e + e~ annihilation, there are only four diagrams at leading order, two of which 
contribute to the fragmentation topologies. On the other hand, the dynam- 
ics underlying the 3 5i-quarkonia production at the Tevatron is much more 
involved. The number of Feynman diagrams is significantly larger (42 vs. 4) 
and only a few of them give rise to fragmentation topologies. As a result, the 
fragmentation approximation starts to be relevant at a much higher regime in 
P T . A similar situation holds for the process 77 — > J /ipcc [32] and also for the 
B* hadroproduction, where the fragmentation approximation is not accurate 
in the Pt range explored at the Tevatron [33,23]. We have also shown that the 
J/ip or T produced in association with a heavy-quark pair of the same flavour 
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are unpolarised. Further analysis, including combination with the results of 
Ref . [34] for pp — > Q + X at NLO and Ref . [35] , to compare with the available 
data from the Tevatron, is ongoing. 

In conclusion, we look forward to the measurement of the fraction of events in 
the J / if) sample at the Tevatron, with at least one charmed meson in the final 
state. Such a measurement could also provide further insight to the mechanism 
responsible for inclusive heavy-quarkonium production at hadron colliders. 
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